Frontier in Astroparticle Physics and Cosmology 



1 



Cosmological Origin of Small-Scale Clumps 
and DM Annihilation Signal 

Veniamin BEREZINSKY 

Laboratori Nazionali del Gran Sasso, Istituto Nazionale di Fisica Nucleare 
O : 1-67010 Assergi (AQ), Italy 

^ : Vyacheslav DOKUCHAEV 

O , Institute for Nuclear Research of the Russian Academy of Sciences, 

60th Anniversary of October Prospect la, Moscow 117312, Russia 
Yury EROSHENKO 

Institute for Nuclear Research of the Russian Academy of Sciences, 
60th Anniversary of October Prospect 7a, Moscow 117312, Russia 



O 



> 
in 
o 
m 

CN ■ Abstract 

^ ! We study the cosmological origin of small-scale DM clumps in the hier- 

r-| ! archical scenario with the most conservative assumption of adiabatic Gaussian 

^! fluctuations. The mass spectrum of small-scale clumps with M < IO^Mq is 

O \ calculated with tidal destruction of the clumps taken into account within the hi- 

\ erarchical model of clump structure. Only 0.1 — 0.5% of small clumps survive the 

stage of tidal destruction in each logarithmic mass interval A In M ~ 1 . The mass 
distribution of clumps has a cutoff at M^am due to diffusion of DM particles out of 
^ ! a fluctuation and free streaming at later stage. Mmin is a model dependent quan- 

5^ \ tity. In the case the neutralino DM, considered as a pure bino, Mmin ~ 1O~^M0. 

The evolution of density profile in a DM clump does not result in the singularity 
because of formation of the core under influence of tidal interaction. The radius 
of the core is Rc ~ O.li?, where R is radius of the clump. The applications for 
annihilation of DM particles in the Galactic halo are studied. The number density 
of clumps as a function of their mass, radius and distance to the Galactic center 
is presented. The enhancement of annihilation signal due to dumpiness, valid 
for arbitrary DM particles, is calculated. In spite of small survival probability, 
the global annihilation signal in most cases is dominated by clumps, with major 
contribution given by small clumps. The enhancement due to large clumps with 
M > lO^M© is very small. 
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1. Introduction 

The gravitationally bound structures in tlie universe are developed from 
primordial density fluctuations 6{x,t) = Sp/p. They are produced at inflation 
from quantum fluctuations. The predicted power spectrum of these fluctuations 
has a nearly universal form P{k) oc with Up ~ 1. At radiation-dominated 
epoch the fluctuations grow logarithmically slowly. After transition aX t = tcq to 
the matter-dominated epoch, the fluctuations grow as 5 oc t'^^^. The gravitation- 
ally bound objects are formed and detached from cosmological expansion when 
fluctuations enter the non-linear stage S > 1. The non-linear stage of fluctuation 
growth has been studied both by analytic calculations P and |2] and in numerical 
simulations [21111111 for Large Scale Structure (LSS). The density proflle in the 
inner part of these objects is given by p(r) oc r~^, with with 1.7—1.9 in ana- 
lytic calculations /5 = 1 in simulations of NFW |2] and (3 = 1.5 in simulations 
of Moore et al. |1] and Jing and Suto [HI. In this work we apply this approach to 
the smallest DM objects in the universe, which we shall call clumps. The clumps, 
being the smallest structures, are produced flrst in the universe, and it makes 
difference of our consideration from LSS formation. The theoretical observation 
of this work is the importance of tidal interaction in the process of DM clump 
formation: the central nonsingular core is formed in the clumps and large fraction 
of clumps are tidally disrupted. 

We use in the calculations the hierarchical model, in which due to merging 
of objects a small clump is hosted by the bigger one, the latter is submerged to 
more bigger etc. We use the standard cosmology with WMAP parameters. The 
primordial spectrum index is Up = 0.99 ± 0.04 (WMAP) or Up = 0.93 ± 0.03 
(WMAP+2dF+Lya). 

2. Tidal Destruction of Clumps in Hierarchical Model 

The destruction of clumps by the tidal interaction occurs at the epoch of 
their hierarchical formation, long before the formation of galaxies. This interac- 
tion arises when two clumps pass near each other and when a clump moves in 
the external gravitational fleld of the bigger host to which this clump belongs. 
In both clump is exited by the external gravitational fleld, i. e. its con- 

stituent particles obtain additional velocities in the c. m. system. The clump 
is destroyed if its internal energy increase AE exceeds the corresponding total 
energy \E\ ~ GM^/2R. In ^ we have calculated the rate of excitation energy 
production by both aforementioned processes. The dominating process is given 
by tidal interaction in the gravitational fleld of the host clumps, with the main 
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contribution from the smallest host clump. We use the Press-Schechter formal- 
ism jHl for hierarchical clustering. A small-scale clump during its life can be a 
constituent part of many host clumps of successively larger masses. After tidal 
disruption of the lightest host, a small clump becomes a constituent part of the 
larger host etc. Transition of a small clump from one host to another occurs 
nearly continuously in time up to the formation of a big enough host, where tidal 
destruction becomes inefficient. 

The fraction of mass in the form of clumps which escape the tidal destruc- 
tion in each logarithmic mass interval A In M ~ 1 is found as 

ei„t^0.01(n + 3). (1) 

In other words the mass function of clumps is C,intdM/M. Since n is close to 
—3, only a small fraction of clumps about 0.1 — 0.5% survive the stage of tidal 
destruction. However, this fraction is enough to dominate the total annihilation 
rate in the Galactic halo. 



3. Core of Dark Matter Clump 

We use the following parameterization of the density profile in a clump: 



Pmt[r) 



Pc[^) , Rc<r<R; (2) 
0, ^ r>R. 



In 12] the relative core radius of the clump is estimated as Xc = Rd -R — -C 1 
from consideration of the perturbation of the velocity field due to damped mode 
of the cosmological density perturbations. Here (5eq is an initial density fiuctuation 
value at the end of radiation dominated epoch. In jTUj the core is considered to be 
produced for spherically symmetric clump by inverse fiow caused by annihilation 
of DM particles. We show [2] that these phenomena are not the main effects and 
that much stronger disturbance of the velocity field in the central part of clumps 
is produced by tidal forces. The tidal forces infiuence the nearly radial motion of 
DM particles at the time of clump formation. As a result these particles obtain 
some angular momentum which prevents the formation of singularity. Once the 
core is produced it is not destroyed in the evolution followed. The core formation 
proceeds mainly near the time of the clump maximal expansion ts- At this moment 
the clump decouples from an expansion of the universe and contracts in the non- 
linear regime. Soon after this period a clump enters the hierarchical stage of 
evolution, when the tidal forces can destroy it, but surviving clumps retain their 
cores. 
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The calculations proceed in the following way (see [7] for details). The 
background gravitational field (including that of the host clumps) is expanded 
in series in respect to the distance from the point with maximum density in a 
fluctuation. The motion of a DM particle in this field is studied. The spherically 
symmetric term of the expansion causes the radial motion of a particle in the 
oscillation regime. Spherically non-symmetric term describes the tidal interaction. 
It results in deflection of a particle trajectory from a center (point with maximum 
density). The average (over statistical ensemble) deflection gives the radius of the 
core Rc- After statistical averaging, Rc is expressed through the amplitude of the 
fluctuation 5cq and the variance (Xcq (or v = 5cq/<7oq) as 

x, = RjR^0.3u-^fi6,^), (3) 

where the function /(5oq) ~ 1 is given in Ref The fluctuations with u ~ 
0.5 — 0.6 have Xc ~ 1, i.e. they are practically destroyed by tidal interactions. 
Most of galactic clumps are formed from ~ 1 peaks, but the main contribution 
to the annihilation signal is given by the clumps with u ~ 2.5 for which Xc — 0.05. 

4. Clumps of Minimal Mass 

The mass spectrum of clumps has a low- mass cutoff at M = M^i^, which 
value is determined by a leakage of DM particles from the overdense fluctuations 
in the early universe. CDM particles at high temperature T > Tj ~ 0.05m^ are in 
the thermodynamical (chemical) equilibrium with cosmic plasma. After freezing 
att > tf and T < Tf, the DM particles remain for some time in kinetic equilibrium 
with plasma, when the temperature of CDM particles is equal to temperature 
of plasma T. At this stage the CDM particles are not perfectly coupled to the 
cosmic plasma. Collisions between a CDM particle and fast particles of ambient 
plasma result in exchange of momenta and a CDM particle diffuses in the space. 
Due to diffusion the DM particles leak from the small-scale fluctuations and thus 
their distribution obtain a cutoff at the minimal mass Md. 

The DM particles get out of the kinetic equilibrium when the energy re- 
laxation time for DM particles Trei becomes larger than the Hubble time H~^{t). 
This conditions determines the time of kinetic decoupling td. At t > td the CDM 
matter particles are moving in the free streaming regime and all fluctuations on 
the scale of free-streaming length A/s and smaller are washed away. The cor- 
responding minimal mass Mfg = (47r/3)p^(to)A|,, is much larger than M^, and 
therefore Mmin = Mfg. In j7] we have performed the calculations using two meth- 
ods: the transparent physical method, based on the description of diffusion and 
free streaming, and more formal method based on solution of kinetic equation 
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for DM particles starting from the period of chemical equilibrium. Both methods 
agree perfectly. The minimal mass in the DM mass distribution is determined by 
the process of free-streaming. For the case of neutralino (bino) as DM particle this 
minimal mass equals to Mmin = 1.5 x 10~^ Mq for neutralino mass = 100 GeV 
and the mass of selectron and sneutrino M = 1 TeV. Our calculations agree rea- 
sonably well with that of [TTj, while from [T2] coincides with our value for 

5. Annihilation Signal Due to Small Clumps 

There is distribution of clumps in the Galactic halo over three parameters, 
mass M, radius R, and distance from the Galactic Center /: nci{M,R,l). This 
distribution also depends on the parameters which describe the internal structure 
of the clumps, P and Xc = Xc{M,R), from Eq. (0). With the number density of 
clumps in the halo written as dN^i = nci{l, M, R)d^ldMdR, the observed signal 
at the position of the Earth from DM particle annihilation in the clumps is given 
by quantity 

/d = — I^CsinC J J dM J dR 

A/jnin Rmin 

X nci(/(C,r),M,i?)iVd(M,/2), (4) 

where r is distance from the Sun (Earth) to a clump and ( is angle between the 
line of observation and the direction to the Galactic center, rmax is the distance 
from the Sun to the halo's outer border and Nci{M, R) is annihilation rate in the 
single clump of mass M and radius R. 

Additional annihilation signal is given by unclumpy DM in the halo with 
homogeneous {i.e. smoothly spread) density PdmIO; where / is a distance to the 
Galactic Center. 

(C) 

/hon. = / rfCsinC / drplMC,r))/ml. (5) 



The enhancement t] of the signal due to a presence of clumps is given by 

V = — f (6) 

-'horn 

This quantity describes the global enhancement of the annihilation signal ob- 
served at the Earth (e. g. the flux of radio, gamma, and neutrino radiations) as 
compared with usual calculations from annihilation of unclumpy DM. We assume 
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Fig. 1. The global enhancement 7] of the annihilation signal from the Galactic halo 
as a function of the minimal clump mass M^[^, for clump density profile with index 
P = 1.5 and for different indices Up of primeval perturbation spectrum. The curves 
are marked by the values of Up. 



that space density of clumps in the halo, nci(/) is proportional to the unclumpy 
DM density, Pdm(0- ^ci(0 = ^Pdm(0/^^ with ^ ^ 1. This assumption holds 
with a good accuracy for the small-scale clumps. The signal from small clumps 
is determined mainly by clumps of the minimal mass. In calculations [7] we used 
different density profiles in the clumps, the distribution of DM clumps over their 
masses M and radii R, and the distribution of clumps in the galactic halo. The 
enhancement depends on the nature of DM particle only through M^^^. For 
(3 = 1.5 enhancement is given in Fig. El Details of calculations and the plots for 
other values of (3 can be found in [7j. The enhancements rj for rip = 1 or less 
is not large: typically it is 2-5 for Mmin ~ 10~^Mq. For example, rj = 5 for 
Up = 1.0 and Mmin = 2 ■ 10~^Mq. It strongly increases for smaller M^in and 
larger Up. For example, for Up = 1.1 and Mmin = 2 ■ 10~^Mq, enhancement is 
?7 = 130. Our approach is based on the hierarchical clustering model in which 
smaller mass objects are formed earlier than the larger ones, i. e. (Teq(M) dimin- 
ishes with the growing of M. This condition is satisfied for objects with mass 
M > M^in — 2 ■ 1O~^M0 only if the primordial power spectrum has the power 
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index Hp > 0.84. The enhancement of the annihilation signal is absent, e.g ?7 — 1, 
for Up < 0.9. 

6. Annihilation Signal Due to Big Clumps 

Numerical simulations reveal in the galactic halos the big clumps with 
masses 10*^ — IO^^Mq. At t ~ t^q these clumps are characterized by an effective 
power spectrum P{k) oc /c" with n ^ —2 (in contrast to n ~ —3 for small clumps), 
and thus the survival probability given by Eq. is larger for the big clumps. 
Indeed, the effective power spectrum index 



n = —3 



1 + 2 



glnaeq(M) 
dlnM 



(7) 



tends to n = Up — 4 = —3 for small M and to n ^ —2 for M in the interval 
10^ - IO^Mq. On the other hand the mean internal number density of DM 
particles in the big clumps is much smaller in comparison with that in the small 
clumps, and it compensates the first effect. 

The number density n{M) of the big clumps with M > IO^Mq in the 
numerical simulations found |1] as n{M)dM oc dM/M"' with 7 ^ 1.9 — 2.0 and 
with a mass fraction of clumps e ~ 0.1 — 0.2. Observations of halo lensing give 
smaller values e ~ 0.06 — 0.07. It is interesting to note that the mass function of 
clumps, obtained from Eq. (Q), is close (including the normalization coefficient) to 
that obtained in the numerical simulations for big clumps with mass M > IO^Mq. 
Strictly speaking our calculations are not valid for big clumps, because of their 
destruction in the halo up to the present epoch to and accretion of new clumps 
into the halo. Nevertheless, for the small interval of masses, where the power-law 
spectrum can be used as a rather good approximation, our approach appears to 
be roughly valid. 

Calculations of the enhancement factor 77 from simplified Eqs. and ^ 
are performed by using e = 0.1, the number density distribution of clumps in the 
Galaxy nc\{l) oc Pdm(O) internal density distribution of the DM particles in 
clumps 

p{r) = p,{r/a)-^{l + r/a)-'' (8) 

valid down to the core radius Re- The core is defined as p(r) = Pc = const at 
r < Re- The NEW profile has (3 = 1 and k = 2, while the Moore et al. profile 
has (3 = K = 1.5. 

The other important parameter is clump radius -R, which determines the 
average density p at the epoch of DM virialization in the clump. This quantity is 
calculated from the value of overdensity 5 at the linear stage of clump formation. 
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The values of S have Gaussian distribution and the normahzation of fluctuation 
spectrum was performed as usual to the value of r.m.s. fluctuation at the 8 Mpc 
scale cTg — 1. This approach corresponds to the picture that the big clumps in 
the Galactic halo are similar to the small protogalaxies (galactic building blocks), 
which escape the tidal destruction when capturing by the Galaxy. The tidal 
stripping of the outer parts of the clumps change their structure. Nevertheless 
this process is not important for clumps with mass M <^ IO^^Mq. For minimal 
and maximal masses of the big clump we use Mmin = lO^M© and Mmax = lO^^M©. 

As a result of calculations, identical to those in Section 5, we have found 
the enhancement factor rj = 1.02 for the NFW profile |SI with 7 = 1.9, e = 0.1, 
Rc = a, R/a = 5 (the case of existence of the core), and rj = 1.14 for the case of 
absence of the core Rc = 0. 

For the Moore et al. profile jlj, ?7 = 1.06 for e = 0.1, Rc/a = 0.5, R/a = 5 
and f] = 1.16 for a smaller core Rc/a = 0.2. 

Diminishing of M^i^ increases the enhancement weakly, approximately as 
r){M^in) oc M-°,35. 

We conclude that enhancement of the annihilation signal due to the big 
clumps is small. 

We are grateful to Ben Moore for providing us with the data which allow 
us to normalize the density distribution of DM in the clumps. 
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